Evidence for Instanton-Induced Dynamics, from Lattice QCD 
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We perform a study of the non-perturbative dynamics of the hght-quark sector of QCD, based 
on some recent resuhs of lattice simulations with chiral fermions. We analyze some correlators 
that are designed to probe the Dirac structure of the quark-quark interaction at different scales. 
We show that, in the non-perturbative regime, such an interaction contains very large scalar and 
pseudo-scalar components. We observe quantitative agreement between lattice QCD results and 
the predictions of the Instanton Liquid Model (ILM). Moreover, we study how the the quark-quark 
interaction is modified, when quark loops are suppressed. We observe a dramatic effect related to 
the loss of unitarity, which is naturally explained in the instanton picture. Such an effect cannot 
be explained in a Dyson-Schwinger Equations (DSE) approach, if one assumes a vector quark-gluon 
coupling. Therefore, from the present study it emerges that DSE models with such an ansatz for 
the vertex function are not consistent with QCD. 



The physics of the light-quark sector of QCD is 
strongly influenced by the non-perturbative structure 
of the vacuum, and in particular by the spontaneous 
breaking of chiral symmetry (SCSB). Hence, identifying 
the microscopic mechanism responsible for such a phe- 
nomenon represents a fundamental step toward the com- 
prehension of the strong interaction. Unfortunately, this 
dynamics resides in the non-perturbative sector of QCD 
and has not been completely understood from first prin- 
ciples. 

The typical energy scale of phenomena related to the 
breaking of chiral symmetry is 47r/7r ~ 1.2 GeV , is con- 
siderably larger than the confinement scale, h^QCjj ~ 
lfm~^ ~ 0.2 GeV. From a theoretical perspective, such a 
separation of scales is crucial, because it justifies attempt- 
ing model descriptions of the non-perturbative physics of 
SCSB, without needing to simultaneously take into ac- 
count the dynamics of quark confinement. The common 
feature in all these semi-phenomenological approaches is 
a strong attraction in the flavor-singlet 0+ channel, lead- 
ing to a quark condensate. On the other hand, some of 
the models which have been proposed in the literature 
rely on drastically different microscopic descriptions of 
the non-perturbative quark-quark interaction. 

Historically, the first attempt to explain the breaking 
of chiral symmetry predates QCD and was developed 
in the Nambu-Jona Lasinio model 0], where a chirally- 
symmetric effective Lagrangian, characterized by a scalar 
and pseudo-scalar four-fermion interaction, was postu- 
lated. Later on, the same structure was recovered in 
the context of the ILM|^. The latter approach has the 
advantage of being formulated in terms of quark and 
gluon degrees of freedom and to be motivated from QCD 
through a semi-classical argument. Moreover, it explains 
in a very natural way the structure of the spectrum of 
lowest-lying eigenvalues of the Dirac operator. 

An alternate model description of the non-perturbative 




FIG. 1: The chirality-flip ratio, i?^^(r), in lattice and in two 
phenomenological models. Squares are lattice points of ,8'j. 
Stars are RILM points obtained numerically from an ensem- 
ble of 100 instantons of 1/3 fm size in a 5.3 x 2.65^ fm* box. 
The solid line is the contribution of a single-instanton, calcu- 
lated analytically in |^. The dashed curve was obtained from 
two free "constituent" quarks with a mass of 400 MeV. Such 
a curve qualitatively resembles the prediction of a model in 
which chiral symmetry is broken through a vector coupling 
(like in present DSE approaches). 



sector of QCD which encodes the physics of SCSB has 
been developed in the context of Dyson-Schwinger Equa- 
tions (DSE). In such an approach, one parametrizes the 
low-energy behavior of QCD through an ansatz of the 
infrared structure of the quark-gluon vertex and of the 
propagators 4]. DSE are then solved numerically, in a 
given truncation scheme. 

Although both DSE on the one hand, and the ILM 
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on the other hand, give comparable phenomenology in 
the light hadron sector, they rely on drastically different 
microscopic pictures of the non-perturbative interaction, 
at the 1 GeV scale. Most applications of the DSE devel- 
oped so far assume a simple vector ansatz^ for the quark- 
gluon vertex function, 1x7^. In other words, the non- 
perturbative dynamics is mediated by the exchange of 
one (albeit non-perturbative) gluon at the time. 

On the other hand, in the instanton picture, the non- 
perturbative dynamics is dominated by the 't Hooft in- 
teraction. Through standard bosonization of the 't Hooft 
vertex, such an instanton-induced interaction can be 
thought as being mediated by fields carrying the quan- 
tum numbers of scalar and pseudo-scalar bosons. 

The goal of this Letter is to identify which one of these 
two alternate microscopic pictures is closer to QCD. We 
shall present evidence for the existence of a large scalar 
and pseudo-scalar component of the non-perturbative 
quark-quark low-energy interaction. This evidence is 
built using some recent results from lattice simulations 
with chiral fermionsQ. On the one hand, these results 
agree on a quantitative level with the predictions of the 
ILM. On the other hand, they rule-out any picture in 
which the non-perturbative quark-quark interaction is as- 
sumed to have a vector structure, like in present DSE 
models. 

Our analysis is based on the study of the flavor Non- 
Singlet (NS) chirality-flip ratio, introduced in [^: 



n^(T) -n^(r) 



(1) 



where (r) and Hg (r) are pseudo-scalar and scalar NS 
two-point correlators^ related to the currents Jt:{t) := 
u{T)i^c,d{T) and J5(t) := u{T)dlT). If the propa- 
gation is chosen along the (Euclidean) time direction, 
A%ip(non-fUp)^'^) represents the probability amplitude 
for s.\qq) pair with iso-spin 1 to be found after a time 
interval r in a state in which the chirality of the quark 
and anti-quark is interchanged (not interchanged). No- 
tice that the ratio R^^{t) must vanish as t ^ (no chi- 
rality flips), and must approach 1 as r — > cxd (infinitely 
many chirality Hips). 

In (6j| it was shown that the correlator is a par- 
ticularly useful theoretical tool for studying the non- 
perturbative dynamics of the light quark sector of QCD. 
In fact, R^^{t) receives no leading perturbative contri- 
bution and probes directly the chirality-mixing interac- 
tion. A spectral analysis of R^^{t) indicated that the 
such an interaction is mediated by topological fields. In 



particular, it was found that the rate of chirality flips 
in a quark-antiquark system is proportional to the mass 
of the 77' meson. Moreover, below we shall see that 
R^^{t) is very sensitive to the Dirac structure of the 
non-perturbative quark-quark interaction. 

The NS scalar and pseudo-scalar two-point func- 
tions appearing in |^ have been first calculated in the 
quenched approximation by the MIT group ^ with Wil- 
son fermions and more recentlyby one of the authors, us- 
ing chiral (overlap) fermions [3. The curve for R^^{t) 
obtained from the result of latter calculation'^ are the 
square points plotted in Fig. ^ 

As expected, the lattice data interpolate between and 
1. Notice that the curve has a maximum at about 0.7 fm, 
where the ratio is considerably larger than one. This im- 
plies that, after few fractions of a fermi, the quarks are 
more likely to be found in a configuration in which their 
chiralities is flipped, than to be found in their initial con- 
figuration. Below we shall see that the presence of such a 
maximum is a signature of a chirality mixing component 
of the quark-quark effective interaction vertex. 

We recall that these lattice results have been obtained 
in the quenched approximation. It is important to ask 
what differences should be expected in full QCD. Using 
general QCD inequalities j^, it is immediate to show that 
i?^-^(T) > 1 if and only if Usir) < 0. The negativity of 
such a two-point function represents a severe failure of 
the quenched approximation which appears only at suffi- 
ciently small values of the quark mass. (In the large mass 
limit, the quenched approximation becomes exact.) It is 
a reflection of the fact that, in the quenched approxima- 
tion, the unitarity of the theory is lost. 

In terms of chirality flipping amplitudes, we see that 
the ^s(t) > constraint implies that quarks must 
never be more likely to be found in the flipped chiral- 
ity conflguration than that in the original conflguration, 
^/Hp(t) < Anon-fiipiT). Heucc, wc Can conclude that 
the fermionic determinant suppresses some chirality flip- 
ping events, which are otherwise allowed in the quenched 
approximation. Indeed, the correlators appearing in lfT|l 
have recently been evaluated in full QCD, with Wil- 
son fermions ^^l- It was observed that the condition 
^s{t) > (or, equivalently, R^^{t) < 1) is restored 
in going from quenched to full QCD. Such a dramatic 
qualitative difference between quenched and full QCD 
calculations of ^ can be used to test phenomenological 
descriptions of the non-perturbative dynamics. Indeed 
any realistic model must reproduce a dramatic enhance- 
ment of the chirality flipping amplitude, when quark loops 
are suppressed. 



^ For an example of a DSE model with a more general ansatz for 
the quark-gluon vertex, see Fischer and AlkoferQ. 

^ Notice that the above correlators are defined in coordinate rep- 
resentation, they are not zero-momentum projections. 



^ We recall that, with overlap fermions, the lattice to continuum 
renormalization factors of the pseudoscalar and scalar correlators 
are equal, and drop out in the ratio. 
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FIG. 2: Suppression of chirality flips due to the topological 
screening induced by fermion-loops in the ILM {Nf = 3). 

Let us now discuss how R'^^{t) looks in the ILM and 
DSE models mentioned above. In both approaches chi- 
ral symmetry is spontaneously broken, albeit through 
two very different microscopic mechanisms. As a conse- 
quence, the quarks acquire dynamically an effective mass, 
which triggers chirality mixing. Due to such a mass gen- 
eration process, in both cases i?^"^(r) interpolates be- 
tween and 1. However, one expects crucial qualita- 
tive differences between the prediction of the ILM and 
that of the DSE models^. Let us first consider the result 
for R^^{t) obtained in the Random Instanton Liquid 
Model^. In such a version of the ILM quarks are assumed 
to propagate in a vacuum populated by an ensemble of 
randomly distributed instanton and anti-instantons with 
a density of n = 1 fm~^ and size p = 1/3 fm. It can been 
shown that the RILM accounts for the 't Hooft inter- 
action to all orders, but neglects quark loops (quenched 
approximation) . 

The RILM prediction for R'^^ir) is presented in Fig.Q 
The agreement with the lattice results is impressive^. It 
is quite remarkable that not only does the RILM curve 
displays a maximum in R^^{t), but also its position 
agrees quantitatively with the lattice results. 

The presence of a maximum in R^^{t), and the sub- 
sequent fall-off towards 1 have a very simple explanation 
in the RILM: if quarks propagate in the vacuum for a 
time comparable with the typical distance between two 
neighbor instantons (i.e. two consecutive 't Hooft inter- 
actions), they have a large probability of crossing the 
field of the closest instanton. If so happens, they must 
necessarily flip their chirality, due to the scalar and pseu- 
doscalar structure of the 't Hooft vertex. So, after some 



More generally, we expect differences between the prediction of 
any model with a vector structure of the vertex on the one hand, 
and any model with a scalar and pseudo-scalar coupling on the 
other hand. 

^ In the following discussion one can regard the RILM as a proto- 
type of a model with a scalar and pseudo-scalar non-local inter- 
action between quasi-zero modes. The strength of the coupling 
is tuned by the instanton density, while the delocalization of the 
quark-quark interaction is provided by the instanton size. 

® The RILM points presented in Fig.lTlhave been obtained using a 
quark mass of about 30 MeV, that is the same as the same bare 
mass used in one of the author's lattice simulations. 



time, the quarks are most likely to be found in the config- 
uration in which their chirality is flipped. On the other 
hand, if one waits for a time much longer than 1 fm, then 
the quarks will "bump" into many other pseudo-particles, 
experiencing several more chirality flips. Eventually, ei- 
ther chirality configurations will become equally probable 
and R^^{t) will approach 1. 

The position of the maximum in R^^{t) carries infor- 
mation about the interplay between one-body and many- 
body effects generating chirality mixing. It is interest- 
ing to compare the above numerical RILM results with 
the single-instanton contribution (solid line in Fig. 
which was derived analytically by one of the authors in 
0. From such a comparison, one can see that one-body 
effects, i.e. chirality flips induced at the level of a single 
interaction, dominate the ratio up to Euclidean times of 
the order of 0.5-0.7 fm. The onset of many-body (many- 
instanton) effects is governed by the numerical value of 
the instanton density: the less dilute is the system, the 
earlier many- instanton effects become important. From 
the agreement between ILM and lattice data one may 
argue that the phenomenological value n ~ 1 fm"** is 
indeed realistic. 

The prediction for R^^{t) would be drastically differ- 
ent in the DSE models and in general in all approaches 
based on a vector quark-gluon coupling. In this case, the 
chirality mixing is only due to the dynamical mass gen- 
eration (i.e. by a genuine many-body effect). In fact, 
unlike in the ILM, a single quark-antiquark interaction 
will not interchange the chirality of quarks, because the 
vector Dirac structure of the quark-gluon vertex is chi- 
rality conserving. As a result, even in the quenched ap- 
proximation, quarks are never more likely to be found 
in the flipped chirality state than in the initial chirality 
state, i.e. R^^{t) < 1 for all r. On a qualitative level, 
the DSE prediction for R^^{t) will be similar to that ob- 
tained considering the propagation of a free but massive 
"constituent" quark and anti-quark pair in the vacuum 
(dashed fine in Fig.P). 

From the above analysis we can conclude that the pres- 
ence of a maximum in the lattice results for R^^{t) im- 
plies that the non-perturbative quark-quark interaction 
contains a strong scalar and pseudo-scalar component. 
In other words, the chirality is mixed already at the level 
of a single quark-quark interaction, and not only through 
many-body effects, such as the mass generation induced 
by SCSB. Therefore, these lattice simulations strongly 
support the ILM picture against the DSE models. 

Additional evidence in this direction comes from com- 
paring quenched and unquenched predictions. In the 
DSE model, neglecting quark loops only affects the speed 
of the running of the coupling, but does not gener- 
ate additional chirality flips. On the contrary, we al- 
ready mentioned that a dramatic qualitative difference 
between quenched and unquenched results is expected 
and observed on the lattice. Such a difference is natu- 
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FIG. 3: Suppression of chirality flipping events, due 
to fermion-loop exchange in the ILM. Circles are RILM 
(quenched) results, squares are IILM (unquenched) results. In 
the unquenched model the unitarity requirement, _R^^(r) < 1 
is restored. 



used these data to test model descriptions of the micro- 
scopic quark dynamics. We have found evidence for a 
large scalar and pseudo-scalar component of the effec- 
tive quark-quark vertex. We have argued that this re- 
sult rules out models in which the quarks couple non- 
perturbatively though a purely vector quark-gluon ver- 
tex. On the contrary, we have observed impressive quan- 
titative agreement between lattice and ILM results. In 
addition, the ILM can also explain how quark loops cause 
a drastic suppression of quark chirality flips. 

It a pleasure to thank G. Ripka, T.Schafer, E. Shuryak, 
and W. Weisc for interesting discussions and for critical 
reading of the manuscript. T. D. and P. F. would like to 
thank respectively the Max-Planck-Institute in Munich 
and Physics Department of the Roma-Tre University for 
their hospitality while this paper was written. This work 
was partially supported by the U. S. Department of En- 
ergy, with grant DE-FG03-95ER40894. 



rally explained in the ILM*". If quark loops are allowed, 
then instantons and anti-instantons can interact through 
fermion exchange. Such an interaction generates an at- 
traction between instantons and anti-instantons leading 
to a screening of the topological charg e il lJ and providing 
a realization of the Leutwyler-Smilga relation (topo- 
logical susceptibility vanishing linearly with quark mass, 
at small quark mass). As a result of such a screening, 
quarks crossing the field of an instanton are very likely 
to find, in the immediate vicinity, an anti-instanton which 
restores their initial chirality configuration^ (see Fig.|21l. 

In Fig. 13 we compare the chirality flipping ratio 
i?^"^(T) obtained from a quenched (RILM) and un- 
quenched (IILM^) version of the ILM. We observe that, 
with the inclusion of the fermionic determinant, the con- 
dition R^^{t) < 1 is restored. We stress that, although 
such a restoration must necessarily take place in QCD, it 
represents a remarkable success of the ILM, which is not 
a unitary field theory. 

In conclusion, we have presented a study of the mi- 
croscopic structure of the non-pcrturbative interaction 
in QCD, based on the analysis of the results of some re- 
cent lattice simulations with chiral fermions. We have 

^ P.F. acknowledges a clarifying discussion with E. Shuryak on this 
point. 

* From this discussion it follows that, in a correlated instanton 
vacuum, the topological screening length coincides with the typ- 
ical time interval f between two consecutive chirality flipping 
interactions. Indeed, in llll such a screening length was esti- 
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